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Abstract. We present the first detection of weak gravitational shear at infrared wavelengths, using observations of 
the lensing cluster Abell 1689, taken with the SofI camera on the ESO-NTT telescope. The imprint of cluster lenses 
on the shapes of the background galaxy population has previously been harnessed at optical wavelengths, and 
this gravitational shear signal enables cluster mass distributions to be probed, independent of whether the matter 
is luminous or dark. At near-infrared wavelengths, the spectrophotometric properties of galaxies facilitate a clean 
selection of background objects for use in the lensing analysis. A finite-field mass reconstruction and application 
of the aperture mass (M ap ) statistic are presented. The probability that the peak of the M ap detection (■§■ ~ 5), 



arises from a chance alignment of background sources is only ~ 4.5 x 10" 
singular isothermal sphere model for the cluster is <tiq 
ratio of ~ 40 Mq / Lq inside a 0.44 Mpc radius. 
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. The velocity dispersion of the best- fit 
1 , and we find a K-band mass-to-light 
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1. Introduction 

A cluster acting as a weak gravitational lens distorts the 
shapes of background galaxies by virtue of its tidal grav- 
itational field. The seminal paper of Kaiser & Squires 
(1993) describes how to use this information to obtain 
a parameter-free reconstruction of a mass distribution. 
Their technique, and variants thereof, has been applied 
to optical observations of a number of clusters (e.g. Clowe 
et al. 2000; Hoekstra et al. 2000). 

Prior to this work, there has been no published re- 
port on the detection of gravitational shear in infrared 
observations. For a given integration time, the attainable 
number density of background galaxies for use in weak 
lensing studies is usually much lower than for optical ob- 
servations. In addition, infrared detectors currently lag be- 
hind optical detectors in their field-of-view making them 
even less efficient for wide-field studies. Nevertheless, with 
developments in detector technology, and new wide-field 
infrared cameras becoming available, this is an avenue to 
be explored. Gray et al. (2000) have considered the mag- 
nification induced depletion effect in near-infrared CIRSI 
observations of the cluster Abell 2219, in order to ob- 



tain best-fit SIS (singular isothermal sphere) and NFW 
(Navarro, Frenk & White 1996) models. 

Here we report on J, H and Ks infrared observations 
in which we detect the gravitational shear signature of the 
cluster Abell 1689 (z = 0.182), manifest in the distorted 
images of background galaxies. We start by outlining the 
observations, and the strategy for obtaining a catalogue of 
background galaxies for use in the weak lensing analysis. 
Then we describe how this catalogue was used to perform 
a mass reconstruction of the cluster. The M ap statistic of 
Schneider (1996) is then calculated, along with a quanti- 
tative measure of its significance. Next we use Ks selected 
cluster members to determine the luminosity of the lens, 
and obtain its mass-to-light ratio by comparison with the 
mass derived from fitting an SIS profile to the weak lensing 
data. We conclude with a brief discussion of the results, 
and the implications for future work. For quantitative esti- 
mates, we assume a cosmology with = 0.3, A = 0.7 and 
Hq = 70 krns" 1 Mpc -1 , and that the background galaxy 
population is at z = 1.0. The conversion between angular 
and linear scales at z = 0.182 is l'.O = 0.18 Mpc. 
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Table 1. The integration time (t), image quality (IQ), 
sensitivity limit (SL) and extinction (A) are shown for 
each of the filters with which A1689 was observed. 



Filter 


t (ks) 


IQ (") 


SL (mag) 


A (mag/airmass) 


J 


14.4 


0.77 


25.59 


0.07 


H 


10.8 


0.85 


24.21 


0.04 


Ks 


15.6 


0.73 


23.57 


0.07 



2. Observations and data reduction 

Our data were obtained on the nights 2000 April 29th and 
30th, using the Sofl camera (Moorwood et al. 1998) on the 
ESO-NTT telescope at Cerro La Silla Observatory, Chile. 
The instrument was used in "Large Field Imaging" mode, 
where the field of view is 4'.9 on a side, and the pixel scale 
is O'.'29/pixel. Abell 1689 was observed in three filters, J, 
H and Ks. Individual integrations lasted between 10 and 
30 s and several of these were co-added to form an indi- 
vidual image. Between images, the telescope was moved 
by up to 2'.0, while always keeping the target within the 
field of view. This offset is somewhat larger than is typi- 
cal in such observations; however, it was necessary given 
the high density of galaxies in the centre of A1689. The 
field of view is effectively increased to ~ 7'.0 (1440 pixels). 
During the observations the seeing was typically ~ C/8 
and conditions were perfectly photometric, except during 
the early hours of the first night. 

The data were reduced in the standard way. From each 
image, the zero-level offset was removed, the flat field cor- 
rection was applied, and then an estimate of the sky from 
other images in the sequence was subtracted. The images 
were combined so that the 4L (signal-to-noise) ratio of the 
final image was maximised and with linear registration 
relative to the first frame in each sequence. 

The photometric calibration was done by observing 
several stars from the list of HST NICMOS photometric 
standards (Persson et al. 1998). Although there was some 
thin cirrus early in the first night, we were able to observe 
a sufficiently large number of standards to derive accurate 
instrumental zero points and atmospheric extinction coef- 
ficients in each of the three filters. The total integration 
times (t), the final image quality (IQ), the sensitivity limit 
(l-cr noise limit in 1 square arcsecond) (SL) and extinc- 
tion coefficients (A) of the combined images in the J, H 
and Ks bands are given in Table 1. The instrumental zero 
points were very stable over the two nights and the scatter 
in the fit of the instrumental magnitude versus airmass is 
less than 0.01 magnitudes in all three filters. 

We did not correct for colour terms between the Sofl 
instrumental system and that used for the HST NICMOS 
standards. Preliminary estimates suggest that these terms 
are all less than 0.02 magnitude. Nor did we correct the 
flat fields for illumination effects which cause the photom- 
etry to vary differentially across the field by up to 0.03 
magnitudes. Thus, the accuracy of the relative and abso- 
lute photometry is around 0.03 magnitudes for both. 



3. Obtaining a catalogue of background objects 

The goal described in this section is the acquisition of a 
catalogue of background object positions and ellipticities, 
that can be used in the weak lensing analysis. 

Our strategy was to start by using SExtractor (Bertin 
& Arnouts 1996) to detect objects on each of the J, H 
and Ks images separately. For all filters, an object detec- 
tion and analysis threshold of at least 3 contiguous pixels 
above la was required. Since the J image contained most 
detections, the H and Ks images were transformed to the 
J reference frame, using IRAF tasks (see Tody 1993 and 
references therein). SExtractor was then run in "double- 
image mode", using the J image for detection (with an 
object detection and analysis threshold of at least 2 con- 
tiguous pixels above 0.7c) and each of the J, H and Ks im- 
ages for measurements. For each object, magnitudes were 
measured using both a 15 pixel (4'.'4) diameter aperture, 
and an outer isophote at 70% of the sky noise (J: 24.07, 
H: 23.11, Ks: 22.55 mag/sq. arcsec). Aperture magnitudes 
are used for colour determination, and isophotal magni- 
tudes are used for all other purposes. The LDAC toolsQ 
were used to obtain a single catalogue of these sources. 

Ellipticities for the objects were determined from the 
second moments of their surface brightness using a modi- 
fied version of Nick Kaiser's IMCATf] software. Bright (but 
unsaturated) stars were selected on the basis of their half- 
light radii and were used to determine the point spread 
function (PSF) as a function of position in the images. A 
bi-cubic polynomial was fit to the stellar ellipticities (for 
25 stars) and used to correct the ellipticities of the galax- 
ies using the technique detailed in Kaiser et al. (1995). 
The smearing of the galaxy ellipticities by the now cir- 
cular PSF was removed using the basic method outlined 
in Luppino and Kaiser (1997), but where the pre-seeing 
shear polarizability tensor was fit as a 5th order poly- 
nomial function of r 9 , the Gaussian smoothing radius at 
which the galaxy achieves maximum significance from the 
background, and the post-seeing ellipticity. The resulting 
shear estimate has been shown to be good to a few per- 
cent accuracy for well behaved PSFs such as those in these 
images (Erben et al. 2001). 

Stellar objects were filtered out, based on their locus 
in a magnitude-radius plot. Next, making use of the small 
dispersion in colour at a particular magnitude for cluster 
galaxies, likely cluster members were excised to exclude 
even the faintest cluster members, where the dispersion in 
the cluster sequence becomes larger. 

Further, only galaxies fainter than mj = 16.6 and hav- 
ing a jr > 5 in J were retained. The limit corresponds 
to mj = 23.2, for the smallest objects of similar size to 
the PSF. Finally, only galaxies redder than the cluster se- 
quence galaxies were kept in the catalogue (they are most 
probably at higher redshifts than the cluster; see for ex- 
ample the colour evolution plots in Poggianti 1997). 



1 Freeware available from ftp://ftp.strw.leidenuniv.nl/pub/ldac/ 

2 http:/ /www. ifa.hawaii.edu/~kaiser/imcat 
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Fig. 1. The colour-magnitude diagram for objects selected 
as background galaxies for use in the lensing analysis 
(open stars). The dashed line at rriKs = 18.5 indicates 
the limit down to which objects selected as cluster galax- 
ies were used in the luminosity estimation (solid trian- 
gles). Other objects detected but not used in the lensing 
or mass-to-light analysis, including stars, fainter cluster 
members, and cluster members outside the radius inside 
which the mass-to-light ratio is determined are marked 
with open squares. Here, mx s and mj are aperture mag- 
nitudes, and niKs(iso) is an isophotal magnitude. Note the 
clear colour-magnitude sequence of cluster galaxies. 

The selection process for the catalogue as outlined 
above is the most conservative of our trials, and the fi- 
nal number of sources used in the weak lensing analysis is 
157, just over 3/arcmin 2 . In Fig.|l], these objects are indi- 
cated (open stars) on a colour-magnitude diagram for the 
observations. We experimented with various cuts, which 
induce only a marginal change in the results of the lens- 
ing analysis described in the next Section. 



4. Weak Lensing Analysis 

For the weak lensing analysis, we used the mean of the 
shear estimates determined from the separate J, H and 
Ks images, weighted by the jj of the object detections. 
The best-fit SIS model to the radial shear profile (cen- 
tred on the brightest cluster galaxy) has a velocity dis- 
persion (j id = lOSOl^kms- 1 , with a significance of 6.1 




Fig. 2. The greyscale shows the SofI J-band image, 7'.0 
on a side, N to the top and E to the left. The solid 
contours show the surface mass density k from a finite- 
field mass reconstruction with a smoothing scale of 0'.54, 
plotted at intervals of 0.05 in k relative to the boundary 
value. A contour plot of M> for M ap with a filter scale of 
n = 1'.95 is plotted as dashed contours, corresponding to 
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relative to the null model (ctid = Okms 1 ). This value 
is consistent with the best-fit SIS to WFI optical data 
(<tid = 998^42 kms" 1 ) presented in King et al. (2001). 
The best-fit NFW model has r 200 = 1-94 Mpc and c = 5.7, 
the significance and quality of fit being the same as the 
SIS. We checked that fitting a radial shear profile to the 
objects bluer than the cluster sequence results in an in- 
significant weak lensing signal. 



4.1. Mass reconstruction 

We use the non-parametric finite-field mass reconstruction 
algorithm detailed in Seitz & Schneider (2001). The recon- 
struction was performed with a smoothing scale of 0'.54. 
The result is shown as solid contours in Fig. normalised 
(necessary because of the mass-sheet degeneracy) so that 
the mean surface mass density on a 0'.097 wide border is 
K — 0.1. Note that the mass peak is consistent with that 
of the luminous infrared image. 

Since the noise properties for non-parametric recon- 
structions are not easily interpreted, multiple reconstruc- 
tions were performed during which the catalogue object 
positions were kept fixed and ellipticities were assigned at 
random (with replacement) from the ellipticity ensemble. 
This enabled us to get a handle on the noise properties, 
and confirms that the central mass concentration is not a 
noise artefact. 
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4.2. Aperture mass statistic 

Next we consider the M ap statistic, which measures 
the surface mass density convolved with a compensated 
Mexican Hat filter function U. Equivalently, the suitably 
weighted ellipticities give a practical estimator of M ap 

M ap ( 1 ?) = ^^e ti ( 1 ?)Q(|0 i - 1 ?|), (1) 

i 

where the sum extends over N galaxies, located at po- 
sitions 9i, inside a (filter) radius r\ centred on position 
i9. The quantity et(i9) is the tangential component of the 
ellipticity with respect to position i9, and Q is a filter func- 
tion related to U. We use the simplest form of the filter 
functions from Schneider et al. (1998). 

An asset of this statistic is that the noise properties 
are well understood, and the significance of a M ap de- 
tection is readily calculated by comparing its value with 
that obtained from the same set of galaxies with random 
position angles. The signal to noise ratio j| for M ap (i9) 
is ^|(i9) where a 2 = (Mi) , {•&). The rms 

VeH N * ' randomisations v ' 

dispersion of M ap is calculated in the absence of lensing, 
since the deviation from this approximation is very small 
in the weak lensing regime ; it gives an upper limit on a . 
The of M ap calculated for a filter scale r\ = 1'.95 and 
for 1 000 randomisations, is shown as dashed contours in 
Fig.||. jj ~ 5 for the peak of the detection; its location 
is consistent with the peak of the non-parametric mass 
reconstruction and with the brightest cluster galaxy on 
the infrared image. Also, note that the M ap detection is 
highly significant over a large area. 

The calculation of in the previous paragraph re- 
lies on the fact that M ap (#) has a Gaussian probability 
distribution. In addition, M ap was calculated around the 
peak of the detection for 4 x 10 7 randomisations of the 
background galaxy orientations, and a higher -j| peak was 
obtained for only 18 of these cases. The probability that 
such a random distribution of galaxies would give a j? in 
excess of the observed peak value is therefore ~ 4.5 x 10~ 7 . 
As another check, the position angles of the background 
galaxies were rotated through 45°, and no highly signifi- 
cant M ap detection was obtained. 

5. Mass-to-light ratio 

Selection of cluster members was performed in the Ks- 
band data, down to isophotal itiks = 18.5 and using their 
tight colour-magnitude and colour-colour loci. Further, we 
restricted selection to objects within 500 pixels (2'.43) of 
the brightest cluster galaxy. Outside of this radius, the sky 
noise increases dramatically due to dithering causing a de- 
crease in exposure time. We indicate the objects selected 
as cluster members on Fig.[|. The apparent magnitudes of 
these galaxies were converted to absolute magnitudes, and 
subsequently to luminosities using the evolutionary and 
k-corrections obtained from the single burst, solar metal- 
licity, Scalo (1986) IMF spectral energy distribution con- 
tained in HyperZ (Bolzonella et al. 2000) (amounting to 
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Fig. 3. The luminosity function of the Ks selected cluster 
galaxies, with the best-fit Schechter function indicated by 
the dashed line. Error bars correspond to Poisson noise in 
each bin. 

—0.53 mag) and the extinction correction from our data 
(0.09 mag). 

The differential K-band luminosity function of the 
cluster members is shown in Fig. ||, with error bars corre- 
sponding to the Poisson variance in each bin, along with 
the best-fit Schechter (1976) luminosity function 

•m*.— pi 

The nonlinear least-squares Marquardt-Levenberg algo- 
rithm was used to obtain the best-fit parameters for the 
faint slope a = —1.01 ± 0.14 and characteristic absolute 
magnitude M* = -24.7 ± 0.40 (i.e. m* = 14.6), with a 
reduced x 2 °f 1-3- Taking into account the difference in 
adopted cosmology, this value of M* is very similar to that 
obtained by Barger et al. (1996), for a sample of distant 
clusters. 

The integrated luminosity (extrapolated to zero lumi- 
nosity) for the best-fit Schechter function is 8.0xl0 12 L© 
and the direct sum of galaxy luminosities (down to L m i n ) 
is 8.8xlO 12 L0. These are consistent, given the errors of 
the fit parameters. We use the direct sum in the mass-to- 
light ratio estimate, since we reach well down the lumi- 
nosity function: the fraction of light observed is approxi- 
mately T(a + 2, ^a)/r(a + 2) ~ 97%. In order to esti- 
mate the error on the luminosity, first consider the error 
arising from the various cuts in the colour-magnitude and 
colour-colour planes, given that each object has a pho- 
tometric uncertainty. We started with the original uncut 
catalogue, and created 100 000 mock catalogues, assigning 
the i th object a flux Fj=F i ( mea )-|- gFi( e rr), where g is a uni- 
form gaussian random variable, and F^mea) , Fj( err ) are the 
measured flux and its error, respectively. The magnitude 
and colour cuts were then applied to these catalogues, and 
the dispersion in the recalculated luminosity determined 
to be only ~ 0.35% of the total. Another very small error 
is the Poissonian error on the total flux, which amounts to 
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~ 0.05%. The dominant error is the 2% error on the lumi- 
nosity because of uncertainty in the Ks zero point. Hence 
for the cluster luminosity we take 8.8 ± 0.18 x 10 12 L Q 
(which is a lower limit because ~ 3% of the luminosity is 
expected to be below the Ks magnitude limit, and light 
outside the isophotal radii of galaxies is not accounted 
for). 

The projected mass of the best-fit SIS to the in- 
frared data within the same 2'.43 radius is M=3.44t°'5i x 
10 14 M Q . This gives an upper- limit K-band mass-to-light 
ratio of 39.1 ± 5.7 M Q /L Q . The best-fit NFW model has 
a mass at this radius which is ~10% higher. 

6. Discussion and conclusions 

We find that it is possible to fit a radial shear profile (the 
best-fit SIS model has am = 1030tgokms _1 ), perform a 
mass reconstruction, and obtain a M ap j| ~ 5 detection 
using only the infrared data on a relatively small field. The 
locations of the peaks of the reconstruction and of the M ap 
detection are consistent with that of the brightest cluster 
galaxy in the infrared image. 

The K-band luminosity function of the cluster galax- 
ies is well described by a Schechter function, with a faint 
end slope a = —1.01 ± 0.14 and characteristic absolute 
magnitude M* = —24.7 ± 0.40. In addition, we find a K- 
band mass-to- light ratio of ~ 40 M /L© in the cosmology 
adopted, in accord with the study of rich Abell clusters 
made by Uson & Boughn (1991). 

The capability to perform routine wide-field infrared 
imaging is just around the corner with projects such as 
the UKIRT WFCAM, the CFHT WIRCAM and VISTA. 
Weak gravitational lcnsing observations of distorted galax- 
ies in the infrared will provide another window through 
which to probe the mass distribution on cluster scales. 
Although the number density of background sources 
achievable is roughly an order of magnitude lower at in- 
frared wavelengths, there are a few points worth noting: 

— Colours of galaxies in the near-infrared are dominated 
by light from old evolved stars, so the well defined 
cluster sequence in the colour-magnitude diagram con- 
tains both ellipticals and (unlike the optical cluster 
sequence) spirals. This makes excising potential clus- 
ter members easier using infrared colours, as anything 
redder than the sequence is likely to be at a higher 
rcdshift, and anything bluer is probably at a lower red- 
shift. The infrared selection of background objects can 
then be used as a consistency check of the optical selec- 
tion, to assess the contamination from cluster members 
and foreground objects. 

— With optical observations, measurement of the cllip- 
ticities of high redshift galaxies can be hampered by 
the presence of blue knots of star formation. 

— The total K-band luminosity of the cluster members 
can be compared with a weak lensing mass, in order to 
obtain a mass-to-light ratio. Further, Broadhurst et al. 
(1992) underlined that the K-band luminosity is a good 



measure of the underlying stellar mass irrespective of 
how it assembled. If a high enough resolution mass map 
became available, we could investigate the correlation 
between any peaks in the total mass and stellar mass 
content. 

Although infrared observations are not preferable to 
optical ones if the goal is simply to do a lensing analy- 
sis, this work highlights the complementary nature of this 
passband. Our next step is to combine the infrared and 
optical information. On a galaxy-by-galaxy basis, the se- 
lection of background galaxies can be examined. In addi- 
tion, we can compare the ellipticities of the lensed galaxies 
determined from each data set. Finally, we will consider 
constraints from strong lensing features to obtain the most 
precise description of the mass distribution of this spec- 
tacular lens. 
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